Impact of Nucleon Mass Shift on the Freeze Out Process 
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The freeze out of a massive nucleon gas through a finite layer with time-like normal is studied. 
The impact of in-medium nucleon mass shift on the freeze out process is investigated. A consider- 
able modification of the thermodynamical variables temperature, flow-velocity, energy density and 
particle density has been found. Due to the nucleon mass shift the freeze out particle distribution 
functions are changed noticeably in comparison with evaluations, which use vacuum nucleon mass. 



PACS numbers: 

I. INTRODUCTION 

High-energy nucleus-nucleus collision experi- 
ments are mainly designed for the search and inves- 
tigation of the predicted new state of matter, the 
Quark-Gluon Plasma (QGP), in which quarks and 
gluons would be set free from the color confinement 
observed in normal nuclear matter. Moreover, 
heavy-ion reactions are expected to exhibit other 
phenomena of Quantum Chromodynamics (QCD) 
in the hot and dense environment of the collision 
region, like in-medium modifications of almost all 
hadrons, or the state of Color Superconductivity 
(CSC). In this respect, the nucleus-nucleus colli- 
sion experiments provide a unique way to test the 
validity of present theoretical approaches and mod- 
els of physics of strongly interacting matter. 

On the other side, a characteristic and inevitable 
problem of collision experiments is that in-medium 
modifications of hadrons and the expected new 
states of matter (e.g. QGP, CSC) disappear by 
the end of the reaction. Accordingly, one can not 
directly measure these properties of the produced 
hot and dense medium. Instead, one has to probe 
the initial stages of the collision indirectly by us- 
ing theoretical models to reproduce the observed 
final particle spectra. A detailed understanding of 



the different stages of a relativistic heavy-ion col- 
lision process becomes therefore very compelling. 
The scheme of a representative relativistic heavy- 
ion collision process looks as follows: 

In the very early stage of nucleus-nucleus colli- 
sions, an extremely hot and dense medium is cre- 
ated in which several hundred or even thousands of 
secondary partons are produced. Due to the high 
partonic density, local (perhaps global) thermal 
equilibrium is reached very rapidly, for instance, at 
RHIC or LHC incident energies within (0.3 — 0.5) 
fm/c for gluons, and (0.5 — 1.0) fm/c for quarks 
P, 0, 0) • It has been proposed that since the 
heavy quark flavor production is dominated by the 
relatively slow gluon-gluon fusion, chemical equi- 
libration of the heavy quark flavors (strangeness, 
charm etc.) might stay incomplete during all the 
collision evolution, so that there is a need to im- 
plement a strangeness suppression factor 7s 
Nevertheless, chemical equilibration of gluons and 
light quark flavors is believed to be reached around 
2 fm/c 0. 

In spite of the nature of the produced medium, 
large pressure gradient perpendicular to the colli- 
sion axes drives the system into rapid expansion 
and to cool down. In heavy-ion collisions, below 
the critical temperature Tc ~ 175 MeV several 
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hundred of hadrons emerge forming a strongly in- 
teracting resonance gas. As the fireball cools down 
further, below the chemical freeze out temperature 
Tch, inelastic collisions cease and hadronic abun- 
dances become fixed. This process is usually called 
the chemical freeze out (cFO). Later on, when the 
hadron gas becomes more dilute, below the ther- 
mal freeze out temperature Tt^, the elastic inter- 
actions cease as well. This stage of the collision 
is usually called kinetic freeze out (kFO). Finally, 
the formed hadrons of the thermal freeze out spec- 
trum propagate freely toward the detectors. Re- 
cently, in 7] both the chemical freeze out temper- 
ature Tch and thermal freeze out temperature Tth 
have been determined for several collision scenarios 
and baryon densities. Nonetheless, the sharp dis- 
tinction between chemical and thermal freeze out 
is an idealization, while in a real collision due to 
the short time scales both processes become mixed 
with each other. Therefore, one sometimes calls it 
freeze out process without further specification be- 
tween chemical and thermal freeze out which im- 
plies Tch — Tth- 

Many kind of approaches have been applied for 
the description of the freeze out of strongl y in ter- 
acting matter. Statistical models [tL la. iTfil Ull 
[T2. 13, 14] can reproduce well the measured par- 
ticle multiplicities in most of the collision experi- 
ments done so far. Kinetic models Ts', as well 
as hydrodynamical approaches 17] have proven to 
be able to describe most of the collective phenom- 
ena like the different flow components in heavy-ion 
reactions. However, despite the success in com- 
parison with experiments, the in-medium modifi- 
cations of the hadrons during the freeze out pro- 
cess have not been taken into account yet. In most 
of the former evaluations the vacuum parameters 
of the particles have been implemented. To the 
best of our knowledge the Refs. 0, seem to 
be the only investigations, where the impact of in- 
medium hadron masses (mesons and baryons) on 
the particle ratios during the chemical freeze out 
has been studied. A systematic study about the 
impact of in-medium hadron masses on the kinetic 
freeze out process has not been performed yet. 

But, implementing the vacuum parameters of 
hadrons for describing the kinetic freeze out pro- 
cess is an approximation which may work or may 
not work, depending on the physical system under 
consideration. For instance, both experiments [20l ] 
and theoretical investigations [2lll22 l. '23*] suggest, 
that pions embedded in a hot and dense medium 
suffer only a small mass change. Accordingly, the 
description of the freeze out process of a purely 
pion gas by means of their vacuum parameters 
seems to be a reliable approximation. On the other 
side, the mass of kaons can be shifted considerably 
in a hot and dense medium 24, 25, 26], so that 
taking into account in-medium modifications for 



the kaon component seems to be compelling. 

In this work, we study a nucleon gas and in- 
vestigate how strong the impact of an in-medium 
mass shift of nucleons on the freeze out profile is. 
We compare the obtained results with calculations 
using a vacuum nucleon mass. 

The paper is organized as follows: The freeze out 
process within a finite time-like layer is considered 
in Sec. 2. The nucleon mass shift and it's imple- 
mentation into the freeze out process are outlined 
in Sec. 3. The results of our study are presented in 
Sec. 4. Finally, in Sec. 5 a summary and outlook 
are given. Further notations and a brief mathe- 
matical remark can be found in the Appendix. 



II. FREEZE OUT PROCESS WITHIN A 
FINITE LAYER 

In this section we are focussing on the last stage 
of the collision, the freeze out process, i.e. we start 
our investigation from the time of collision where 
the expanding and cooling down system reaches a 
temperature T < Tc, where the hadronization of 
the primary parton gas is almost completed. 

The frozen out particles are formed in a layer of 
finite thickness L, bounded by two hyper-surfaces: 
the pre- freeze out hyper-surface with TpicFO — 7c 
and a post-freeze out hyper-surface with Tpost fo — 
Tch — Tth- These surfaces are defined by the 
normal dcTp, which in general can be a space-like 
da-fj^da-'^ < or time-like four vector, da^da^ > 0. 
The diameter L of the layer is of the order of a 
few mean free paths of the particle under consid- 
eration. To get an idea about the scales we recall 
that for nucleons at ground state saturation den- 
sity the mean free path is about 1 fm [2^ |2^ . 

Dynamical models, like hydrodynamical or 
transport models, allow to describe such freeze out 
processes through the layer. In doing so, the hy- 
drodynamical models have certain advantages over 
transport model calculations. An important one is 
that, once the equation of state and initial con- 
ditions of the hadronic matter are specified, the 
space-time evolution of the system is uniquely de- 
termined by the hydrodynamic differential equa- 
tions. Especially, this implies that the impact of 
several equation of state may be investigated in a 
very direct way. Even more, uncertainties or as- 
sumptions made in the underlying kinetic theory 
of the particles under consideration are circum- 
vented. In addition, the use of familiar thermody- 
namical concepts, like temperature, flow velocity, 
pressure and energy density also provide a trans- 
parent physical picture of the evolution. Of course, 
the basis of applicability of hydrodynamics is the 
assumption of local thermal and chemical equilib- 
rium. In the following we will suppose the validity 
of these conditions and will apply the theory of hy- 
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drodynamics for describing the thermal freeze out 
process in a finite space-time layer. 

The theoretical description of the kinetic freeze 
out within a hydrodynamical approach has been 

worked out some years ago 29, sdi m m m. 

Very recently, in [SJI and (35|] the formalism has 
been applied to the case of a finite freeze out layer, 
separately both for space-like and time-like nor- 
mals. While the formalism in [s^. Issjl has been 
developed for the general case of a massive par- 
ticle, the calculations have been performed for a 
massless pion gas. Here, we will make use of the 
outlined formalism of Ref. ,35j for a time-like layer 
for the case of a massive nucleon. In particular, we 
will implement the in-medium mass modification 
of nucleons traveling through the freeze out la yer . 
It is not necessary to repeat the formalism of |35l | 
in detail. Instead, we shall restrict our explana- 
tions on the basic concept and will only give the 
equations relevant for our study. 

Local equilibrium implies that the thermody- 
namical parameters inside the layer become space- 
time dependent, i.e. we have a space-time depen- 
dent temperature T{x), flow velocity v{x), energy 
density e{x) and nucleon density n{x). For evalu- 
ating these functions we need the basic equations 
of hydrodynamics, 

d^N^ix)=0, and d^T'^" (x) ^ , (1) 

where 



^k'^f{x,k) 



is the particle current, and 



— fc'^r/(x,A:) 



(2) 



(3) 



is the energy momentum tensor. Here, x'^ = (t, r) 
is the four-coordinate and fc^ = {Ek , k) is the 
four-momentum of the nucleon. While the first 
relation in ^ is only valid when the total num- 
ber of particles is conserved, the second relation in 
is always satisfied and asserts the energy and 
momentum conservation. The one-particle dis- 
tribution function f{x,k) is an invariant Lorentz 
scalar, and is normalized to the invariant num- 
ber of particles N (in our case the nucleons), i.e. 
TV = / d^r d^k f{x,k). 

While the components of the tensors (|2Jl and (j2Jl 
depend on the Lorentz frame chosen, two Lorentz 
invariant scalars can be obtained, the invariant 
scalar energy density e and invariant scalar par- 
ticle density n: 



e{x) = Ufj,{x)T^''{x)u^{x) 



n{x) = Uf,{x)N''{x) 



(4) 



(5) 



We notice that these invariant scalars have to be 
distinguished from the non-invariant energy den- 
sity e — E/V and particle density h — N/V, where 
E is the non-invariant total energy and V is the 
non-invariant volume of the system. 

While the invariant relations Q and (|SJ) are 
valid in any Lorentz frame, in a concrete evaluation 
one has to specify the frame in which the compo- 
nents of the four-current, energy- momentum ten- 
sor and the (always time-like) four-velocity are 
evaluated. Any Lorentz frame can be defined by a 
Lorentz boost in respect to the local Rest Frame 
of the nucleon Gas, RFG, on which the condi- 
tion u'^-pq{x) — (1, 0, 0, 0) is imposed; obviously, in 
RFG we have e = crfg and n = nptPG- However, 
this condition does not define the RFG uniquely. 
There are, in general, several possibilities to de- 
fine such a rest frame. Here, we will take Eckart's 
definition js^, which is the most appropriate one 
for heavy-ion reactions with high baryon densities. 
According to this definition the local Rest Frame is 
tied to conserved particles, which can be achieved 
by equating the unit vector of the particle four- 
current with the four-velocity of the particle flow. 



u^'ix) 



N^'{x) 



y/N-{x)N,{x) 



(6) 



Accordingly, in RFG there is no particle flow in 
spatial directions. It is straightforward to recog- 
nize, that the Lorentz invariant denominator in 
(El is just the invariant scalar particle density of 
Eq. Q. And, while the components of four- vectors 
u'^ and depend on the Lorentz frame cho- 
sen, the tensor relation ©, which connects these 
frame-dependent components, remains valid in any 
frame. 

From the definitions Q , (jsj and © one obtains 
the following set of three coupled differential equa- 
tions, which, by means of Eckart's definition, are 
valid in any Lorentz frame: 



de{x) = u^{x) dT^^ {x) u^{x) 



dn{x) = u^{x)dN^'{x) , 



(7) 



(8) 



duf'ix) 



n{x) 



gf"" -u^'{x)u''ix) dN^ .(9) 



Altogether, since there are four unknowns in the 
problem under consideration, namely T,v,e,n^ an 
additional constraint is necessary to get a com- 
plete system of equations, which uniquely deter- 
mines these four unknowns. That constraint is 
provided by the Equation of State (EoS) for the 
nucleon gas js^. Ha [s^ , which is assumed to be 
valid in any space-time point of the reaction zone 
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after hadronization, 

e{x) = n{x) 

K 

+ 18 



MN{n{x),T{x)) ~ Eo 

2 



n(x) 

HQ 



1 



(10) 



The term Eq — 16 MeV accounts for the nu- 
clear binding energy among the nucleons, the term 
proportional to the compressibility constant K = 
9 {dp/dn)n=na — 235 MeV accounts for the depen- 
dence of compressibility on density. Since we are 
aiming at investigating the impact of in-medium 
nucleon mass shift on the freeze out process, we 
have already implemented a density and temper- 
ature dependent nucleon pole mass in H10() . The 
given EoS (|10|l is a generalization of EoS for the 
ideal nucleon gas, which is valid in the rest frame, 
i.e. eRFG = nRFG[A^Ar + 3/2rRFG]- There are 
other generalizations for the nucleonic EoS |3a |. 
However, we have checked that in the energy and 
temperature region we are working here, the re- 
sults obtained are insensitive on the specific choice 
of nucleonic EoS. The EoS ifTUI) is used to deter- 
mine the temperature T{x) of the interacting com- 
ponent of the nucleon gas during the freeze out pro- 
cess. Accordingly, the four equations 10, ©, © 
and (|10|l represent a closed set for evaluating the 
four unknowns T, w, e, n of the one-particle system. 

Now we will turn to the explicit evaluation of 
components for the energy-momentum tensor and 
nucleon four-current. In line with Ref. [35| we 
will perform all evaluations in the Lorentz Rest 
Frame of the freeze out Front, RFF, so that in our 
study all tensor components in Eqs. 10-10 can be 
labeled by RFF. The Lorentz frame RFF is defined 
as follows: 

At a given instant in the space-time the expand- 
ing hot and dense hadronic system reaches a cer- 
tain freeze out temperature Tpost fo i where all con- 
stituents of the system are assumed to get frozen 
out, i.e. all hadrons do not interact anymore. In an 
arbitrary but fixed direction Gx = rT/|rT| trans- 
verse to the beam the Rest Frame of the gas RFG 
moves with a velocity vt relative to the freeze out 
front RFF. Then, by means of a Lorentz transfor- 
mation the particle four-velocity in RFF becomes 

sign(vT)|vT| and 
2 The velocity v is called flow- 



it 



RFF = where v 

7 = i/\/T 



V 



velocity and, in general, can be positive, negative 
or even zero. 

Furthermore, as the system expands and cools 
down the number of interacting particles decreases 
up to the post freeze out surface of the finite layer, 
where by definition the density of interacting par- 
ticles vanishes. Accordingly, the thermal freeze out 
process inside the layer can be described by decom- 
posing the particle distribution function into two 
components of the matter, an interacting part fi 



and a non- interacting free part //, thus 

f{x,k) = n{x,k) + ff{x,k) . (11) 

According to Eq. ((TT|l and by means of Q, ®, 
© we have an interacting and a non-interacting 
particle density, 

n,{x) = ^N,,{x)N^{x) , 

nf{x) = ^N,f{x)N-ix) , (12) 

with n = Uf + ni. On the pre- freeze out we assume 
to have thermal equilibrium, i.e. we have a Jiittner 
distribution for fi as starting one-particle distribu- 
tion function, while by definition // is zero on the 
pre-freeze out hyper-surface. The space-time evo- 
lution of the interacting and non-interacting com- 
ponents inside the layer is governed by the follow- 
ing differential equations [35j: 



+ -[/e,(i)-/.] , 
To 



(13) 



I f L \ /k^'da^ 



fr, (14) 



with the time r between collisions. The Jiittner 
distribution is given as ^3 



1 



(15) 



with the chemical potential /z; for the inter- 
acting component it is determined by Eq. (|33|l 
given below. The second term in II 1311 is the re- 
thermahzation term [23, 113, 112, IH, js^] which de- 
scribes how fast the interacting component ap- 
proaches the Jiittner distribution within a relax- 
ation time tq. Here, we will use the immedi- 
ate re-thermalization limit tq — > 0, which implies 
fi ~^ feq faster than tq 0, i.e. local equilibrium 
at all times during the freeze out in the following 
way: 

First, the layer is subdivided into small inter- 
vals. Then we calculate the changes dT^'' and dN'^ 
based on their kinetic definitions and (0, re- 
spectively, with the freeze out distribution fi. At 
the beginning of a time step this is considered to 
be a flux coming from a Jiittner distribution, and 
continues during the length of the whole time step 
according to the kinetic differential equation (|13|l 
(without the re-thermalization term). Then the 
remaining distribution is not of Jiittner type any- 
more. Nevertheless, the loss dT and dN are cal- 
culated based on the initial Jiittner and the es- 
cape probability. When we are at the end of the 
time step of such a small interval of the layer, we 
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have a change in all thermodynamical variables 
T,v,e,n. With tq — + we assume an immediate 
re-thermalization of T'^^ and N^^, i.e. we define a 
new Jiittner distribution with the new values for 
T,v,n at the end of the time step. At the next 
time step we use this new Jiittner distribution to 
calculate the changes of it in the next small time in- 
terval, and so on. Accordingly, the last term in ||7J) 
vanishes, as it can be seen as follows. Since at the 
beginning of a time step we take a Jiittner distribu- 
tion according to the immediate re-thermalization 
limit, the second term of Q is zero (see Appendix). 
Then, during a time step the energy momentum 
tensor of a Jiittner distribution is changed by 
an amount of dT^'^ ^ dt, governed by Eq. 
That means the second term in {TJ is of order 
0{du dt), i.e. of second order in the differentials, 
so that the second term in Q has to be neglected. 
For more details about the relations H13|l and H14|l 
and about the re-thermalization limit we refer the 
interested reader to Refs. ^29.. ,30„ ,33„. .SSj . 



By means of the microscopic definitions jSJl and 
(PI one obtains for the change of the four-current 
and energy momentum tensor the following general 
expressions for the interacting component , 



dNt = dt I ^k'' 



drr - dt I ^kf'k" 

k^ 



dtf, 
dtf, 



(16) 
(17) 



Since we are mainly interested on the freeze out of 
the interacting nucleons we will write down the in- 
teracting component of these tensors and drop the 
index i in the following. The non-interacting com- 
ponents can be deduced from them by changing 
the sign in front. We will write down these expres- 
sions explicitly for the change of dN and dT as 
given in Ref.|35l| for the RFF, which, as previously 
mentioned, are related to the RFG by a Lorentz 
boost: 



dN°{t,v,T,MN,n) 





dt 


dN'^it, 


v,T,MN,n) 




dt 


droo(t. 


v,T,MN,n) 




dt 


dT^'^it, 


v,T,MN,n) 




dt 


dT'=''{t, 


v,T,MN,n) 


dt 



^j^^j{G^{MN,v,T)-Gt{MN,v,T)) , 

1 dN"{t,v,T,MN,n) ^ 1 L n { AaKi{a) ^ 2a^Ka{a) 
V dt 

I L nT 1 



T L — t 4 

1 dT°°{t,v,T,MN,n 
V dt 

1 L nTb'^ 



T L-t 4 \ V 
{G^{Mn,v,T) - G+{Mn,v,T)) 



{{■i + v'^)K2{a)+aKi{a)) 



TL-t 2 V 

1 dT°''{t,v,T,MN,n) 
V dt 

T /dN''{t,v,T,MN,n) 1 dN°{t,v,T, MN,n) 



7 V 



dt 



dt 



T L — t 2 V w 



(18) 
(19) 
(20) 



(21) 



(22) 



Here, a = Mn/T and b = "f a. The functions 
and Kn are defined in the Appendix. The 
set of equations - (jTHll and (UHll - IHl) allow 
us to evaluate the basic thermodynamical function 
T{x), v{x), e{x) and n{x) during the freeze out pro- 
cess for a particle with a constant mass Mjy. How- 
ever, as mentioned in the Introduction we are aim- 
ing at an implementation of in-medium mass shift 
to look for it's impact on the freeze out process. 



Therefore, we will first evaluate the equation with 
the vacuum nucleon pole mass Mjv(O), and after- 
wards replace it by a density and temperature de- 
pendent nucleon pole mass {n, T) . 
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III. NUCLEON MASS SHIFT 

During the freeze out process, the temperature 
and particle densities are presumably close to the 
deconfinement phase transition critical values 0|- 
Therefore, the in-medium values of masses, decay 
widths, coupling constants, and all other physical 
quantities characterizing the particles under con- 
sideration have to be taken into account. In our 
study we examine a purely nucleon gas, and con- 
sider the in-medium mass modification of nucleons 
located in a hot and dense nuclear enviroment. 

We start with a brief reconsideration of the nu- 
cleon mass in vacuum. The nucleon derives it's 
vacuum mass, Mn{0) = 939 MeV, from the quark- 
gluon interaction of it's underlying substructure, 
consisting of valence quarks, sea quarks and glu- 
ons. However, although there has been consider- 
able success in reproducing the vacuum mass of 
nucleons on the basis of their microscopic quar k 
and gluon substructure (lattice evaluations, HJ), 
a rigorous use of fundamental theory of QCD in 
this respect is not yet in reach. Therefore, our un- 
derstanding of the nucleon's mass structure comes 
mostly from models. From a hadronic field theo- 
retical point of view the nucleon mass Mn{0) can 
be defined as the pole mass of the nucleon propa- 
gator in vacuum, 

HNik) = i j d^xe'''^ (0|T#jv(a;)^Ar(0)|0) 

= ^ , (23) 

7^^^" Mn -EN{k) + ie 

where T is the Dirac time-ordering, ^'^v is the nu- 
cleon field operator, '^'n — *I'jv7o> In are the Dirac 
matrices, and 'EN(k) is the nucleon self energy in 

o 

vacuum. The parameter M n is called bare nucleon 
mass, i.e. the mass parameter entering the La- 
grangian which describes the interaction between 
the nucleons and other hadrons (e.g. nucleon-pion 

o 

interaction). In general, the mass parameter Mn 
has to be distinguished from the vacuum pole mass 
of nucleon AfAr(O) 939 MeV, defined by 

Mjv(O) =Mn +RcI]Ar(7^fc^ = MAr(O)) . (24) 

As mentioned, there are several models which al- 
low to calculate the pole mass from a QCD based 
microscopic point of view. Among them is the 
extension of QCD sum rule approach to the 
case of baryons ,43] , which provides an interlook- 
ing between the nucleon pole mass and QCD based 
quantities, so called QCD condensates. Within the 
QCD sum rule approach, the nucleon field opera- 
tor ^fjv in is expressed by an interpolating 
field f)jv [4^ , which is made of up and down quark 
field operators, and which has the quantum num- 
bers of a nucleon (charge, spin, isospin, parity). In 



this line, in the so called loffe formula for the 
nucleon pole mass in vacuum has been obtained, 

MNiO) = (0|gq|0) , (25) 

providing a link between the pole mass and the 
chiral condensate, {0\qq\0) = -(0.250 GeV)^; M ~ 
1.15 GeV is the Borel mass parameter determined 
by stab ility constraint of the nucleon sum rule ap- 
proach [43 |. 

Now we will turn to the in-medium nucleon pole 
mass M]\[{n,T), which is the very characteristics 
which enters the EoS 1)10(1 . |4l |. In general, a nu- 
cleon propagating in a hot and dense hadronic en- 
viroment can be regarded as a quasi particle, de- 
scribed by the in-medium nucleon correlator 

njv(A:,?i,T) 

= i J d*x e'^"" {n[liiN{x)^N{'^)\^) , 

= -0 ^ • (26) 

"ff_ikt^- Mn -Sjv(fc, n, T) + ie 

The nucleon self energy TiNik^njT) in medium de- 
pends on density and temperature of the surround- 
ing hadronic medium inside the freeze out layer; 
the hadronic medium is described by the state jfi). 
In generalization of Eq. I|24|) the in-mcdium nu- 
cleon pole mass is defined by 

Mjv(n,T) 

^Mn +ReY.Nh|,k^' = MN{n,T),n,T) . 

(27) 

From this point of view it becomes obvious that 
the pole mass will be modified in a hot and dense 
hadronic matter, simply due to the fact that the 
self energy of a nucleon in medium will be different 
from a nucleon in vacuum. 

In general, the particle pole mass in (|27|l for a 
nucleon at rest (RFC), embedded in a hot and 
dense hadronic medium, is given by [45I |4^ 

M^v(n,r) 

= MAr(O) + ReI]s(n,T) + ReI]v(?i,r) , 

(28) 

with the attractive scalar part (Re E5 < 0) and 
the repulsive vector part (Re > 0) of nucleon 
self energy in medium. It is a result of several 
theoretical models applied so far, that the individ- 
ual contributions of scalar and vector self energy 
are large, but they are canceled by each other to 
a large extent; typical values at saturation den- 
sity are Re Ec = -400 MeV, Re Ey = +300 
MeV m, H 133. In particular, several theoret- 
ical approaches predict a mass dropping of the nu- 
cleon pole mass in a hadronic enviroment of about 
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MN{no,0)-MN{0) ~ -(80 ±20) MeV at ground- 
state nuclear saturation density uq ~ 0.17 fm^'^ 
and at vanishing temperature. Here, we will take 
the QCD sum rule results for a nuclcon in matter, 
given by 13 

ReE,(.,T) = +Af.(0) f - lV(29) 



ReEy(n,T) 



{0\qq\0) 
--Mat OJ 

3 {Om\0) 



(30) 



where we have accounted for the lowest mass di- 
mension condensates only; gluon condensate and 
higher mass dimension condensates give only small 
corrections due to large cancellations between their 
individual contributions. We recall that the part 
= Mjv(O) + ReEs of ^ resembles the termi- 
nology " effective mass" used in the Walecka model 
[isjl . the Skyrme model |4^, and which has also 
been evaluated by means of a mean field approach 
in [s^l ■ For a more detailed clarifying of the ter- 
minology "effective mass" , of ten used in different 
meaning, we refer to Ref . [SJ , where Af^ is called 
Dirac mass. 

The chiral condensate at finite temperature and 
density, (f2|gq|r2), has been evaluated within the 
Nambu-Jona-Lasinio (NJL) model in Ref. [s^. 
Later, in Ref. |53j | the in- medium chiral conden- 
sate has been evaluated at finite densities and tem- 
peratures by means of a pion-nucleon gas, finding 
a good agreement with the results of Ref. ,52,] . 
Here, the condensates (|29|l and H30|l have to be 
evaluated for a purely nucleon gas to be consistent 
within the whole approach presented. According 
to Eqs. lfTT|l and (fT^ there are two components 
inside the finite layer: an interacting component 
with density nt and a non-interacting component 
with density Uf. For evaluating the condensates 
(1^ and we approximate the interacting com- 
ponent by a Fermi gas with chemical potential Hi 
and temperature T. On the other side, the temper- 
ature for the non-interacting component becomes 
ill-defined. Nonetheless, a relevant physical pa- 
rameter for describing the non-interacting compo- 
nent remains the density n f . Accordingly, the con- 
densates in one-particle approximation are given as 
follows ii: 



{n\qq\n) = {0\qq\0} 



+4: 



d^k 1 

(27r)3 2E, 



■NF{N{k)\qq\N{k)) 



{n\qh\n) 



2Mn{0) 



{N{k)\qq\Nik)) , (31) 
NF{N{k)\qU\N{k)) 



(2^)3 2Ek 

-{N{k)\qU\N{k)) 



2Mjv(0) 



(32) 



where Np = [e(^k-tJ;)/T ^ ^-i jg ^^le Fermi 
distribution, and the nucleon energy is Ek = 
^AfAr(0)2 -|-k2. Note that {0\qh\0) = 0. Here, 
the relativistic normalization (7V(ki)| A^(k2)) = 
2i;fej(27r)35(3)(ki - kz) is used. In Eqs. (jSTjl and 
(|23) the spin (up, down) and isospin (proton, neu- 
tron) degeneracy of nucleon states has been taken 
into account by the factor 4 in front of the mo- 
mentum integrals. The chemical potential for the 
interacting component can be evaluated via 



d^k 



1 



(27r)3 e(-Bfc-p.)/T + 1 



(33) 
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FIG. 1: Effective in-medium nucleon pole mass 
Afjv(n, T) according to Eq. H28|l (for more details see 
main text). 



The condensates in Fermi gas approximation are 
given by p3| 



{N{k)\qq\N{k)) = 



MN{0)aN 



{N{k)\qU\N{k)) = 3A/jv(0) 



(34) 



(35) 



The nucleon sigma term is crjv — 50 MeV |55|, and 
TOq ~ 5 MeV is the averaged current quark mass 
of the up and down quark flavor Insert- 
ing these parameters into (|29l) and H3U|) wc obtain 
Re E5 = -390 MeV and Re Sy = +315 MeV at 
ground state saturation density hq. The Eqs. - 
H35|l summarize our propositions made for obtain- 
ing the in-medium nucleon pole mass MM{n,T) 
which enters the EoS (fTUIl . Fig. Q shows the 
dropping of the in- medium nucleon pole mass. The 
slight increase of the in-medium nucleon pole mass 
with temperature is an artifact of the purely nu- 
cleon gas approximation. That means, an imple- 
mentation of pions in Eqs. (|31|1 and H32() which 
govern the mass relation (|28|l . would cause a tem- 
perature decreasing of these condensates js^l and 
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then of the in-medium nucleon pole mass. Here, in 
a baryon dominated system this artificial increase 
of Mi\[{n,T) with temperature is superposed by 
the much stronger down shift of the pole mass with 
nucleon density. 



IV. RESULTS AND DISCUSSION 

In this section we represent and discuss the re- 
sults of the coupled set of differential equations 
O - (jni in combination with the EoS (jTUJ and 
the in-medium nucleon mass shift relations l|28|) - 
(I35|l . The differential equations have been solved 
with the Runge-Kutta method [fisl IHfll |^ on 
the IBM 1300 cluster at Bergen Center for Com- 
putational Science (BCCS). For all of the cal- 
culations, we have taken Tprc fo = 150 MeV, 
ripreFO = 1.5 no (corresponding to /Xprc fo — 
615 MeV) and f pre fo = 0.5 c as starting values 
on the pre freeze out hypersurface. These values 
are, for instance, in line with typical parameters 
which have been reached within the Alternating- 
Gradient Synchrotron (AGS) at Brookhaven Na- 
tional Laboratory (BNL) in Brookhaven/USA, cf. 
|6lj . Higher baryonic densities can be reached 
within the Schwer-Ionen-Synchrotron (SIS) at 
Gesellschaft fiir Schwerionenforschung (GSI) in 
Darmstadt/Germany, cf. [g^- Note that Tpro fo 
and npre fo are pre freeze out values and, there- 
fore, they are larger than typical post freeze out 
values given, for instance, in Ref. j2i]. 

In Figs. |5] and 131 the time evolution of the pri- 
mary thermodynamical functions through the fi- 
nite freeze out layer arc shown, in terms of the 
proper time t. Note that the densities n = rii + n f 
and e = Ci + Cf are kept constant inside the layer. 
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FIG. 2: LEFT: The temperature of the interacting 
component. RIGHT: The flow velocity parameter v 
of the interacting component. The solid lines are with 
a constant nucleon mass Mn{0) = 939 MeV, while 
the dashed curves are evaluated with a density and 
temperature dependent nucleon mass Mjv(n,r). 



We find a substantial impact of in-medium mass 
modification on the freeze out process within the 
purely nucleon gas model. Furthermore, the 
Figs. 121 and also elucidate, that the freeze out 
process proceeds faster for all thermodynamical 
quantities T, v, e, n when taking into account the 
mass dropping of nucleons. The physical reason 
for a faster freeze out originates from a smaller en- 
ergy density of the nucleon system due to a smaller 
nucleon mass Mn {n, T) compared to the vacuum 
nucleon mass M7v(0). 

The given functions for T,v,e,n are not directly 
accessible. In experiments the way to study the 
hot and dense hadronic matter produced in heavy- 
ion collisions is to measure the distributions of final 
state particles, which reach the detectors long time 



after their last interaction. Accordingly, as next 
we consider the one-particle freeze out distribution 
function at ky — 0, i.e. ffikx) = ff{kx,ky = 0) 
and consider the impact of the evaluated thermo- 
dynamical functions T,v,e,n on it. The results 
are shown in Fig. ^ for different instants during 
the freeze out process. The function //(fc^) is de- 
termined at thejpoint A [6^ of the freeze out front; 
see also Ref. ,35| for more details. The function 
ff{kx) is obtained by solving the differential equa- 
tion (|14|l . where for fi the Jiittner distribution (|15|l 
is used, but with the parameters T and v as deter- 
mined previously and given in Fig. |21 

The logarithmic scale in Fig. ^ disguised the 
strong modification of these distribution functions. 
For small moments up to k^ < 1 GeV, at the very 
beginning of the freeze out process at t = 0.1 r 
there is a change of Jf{kx), which remains up to 
the end of the freeze out process aX t = 9.0 r. A 
contour plot of the freeze out particle distribution 
functions ff{kx, ky) over their transversal and Ion- 
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FIG. 3: LEFT: Nucleon energy density of the inter- 
acting component. RIGHT: Nucleon particle density 
of the interacting component. The solid lines are with 
a constant nucleon mass Mjv(O) = 939 MeV, while the 
dashed curves are evaluated with a density and tem- 
perature dependent nucleon mass Mjv(n,T). 



gitudinal momenta and ky, respectively, shown 
in Fig. [SI illustrates this statement. We observe a 
remarkable change by a factor ~ 2 for momenta 
kx ) ky < 1 GeV. 

A few remarks are in order about the used start- 
ing values for density and temperature. First, the 
formulas H29|l and (|30|l have, like other theoretical 
approaches, a limited range of validity in respect to 
the density; n < 1.5 no . And second, according to 
Eq. (|13|l the re-thermalization is assumed within a 
time step dt. Numerical accuracy for solving the 
set of differential equations Q) - © requires suffi- 
ciently small time inter vails dt. However, a smaller 
starting temperature Tpio fo implies a longer re- 
thermalization time To < c?i, so that Tpro po cannot 
be taken arbitrary small. In addition, these two 
boundaries have to be adjusted to be in a region of 
the QCD phase diagram where we are inside the 
hadronization region and above the kinetic freeze 
out. The parameter choice of the starting values. 



FIG. 4: LEFT: Freeze out distribution function 
ff{kx,ky = 0) at different instants t, evaluated with a 
constant nucleon mass Mn{0) = 939 MeV. RIGHT: 
Freeze out distribution function ff{kx,ky — 0) at dif- 
ferent instants t, evaluated with the density and tem- 
perature dependent nucleon mass MMiji, T). The lines 
are as in the left panel. The Freeze out distribution 
function has increased by a factor of ~ 2, when den- 
sity and temperatur dependent nucleon masses were 
taken into account. 



n-ptoFO = 1-5 no, Tpro FO = 150 MeV, are an op- 
timal compromise for these borderlines. Within 
the approach presented we have a common way 
to model the kinetic freeze out process, and which 
is capable to implement the nucleon mass shift by 
means of a purely nucleon gas model. Nevertheless, 
one has also to be aware that the pion-nucleon ratio 
becomes small only for high enough nucleon den- 
sities n = (1.5 — 2)no and moderate temperatures 
T ~ 100 MeV, e.g. Our starting values for 

density and temperature on the pre freeze out hy- 
persurface deviate from these values. Therefore, a 
more sophisticated model requires the implemen- 
tation of pions and maybe even heavier mesons. 
However, due to different freeze out scenarios be- 
tween nucleons and mesons, cf. |6^ . such a proce- 
dure would require the use of a two-fluid or even 
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FIG. 5: Freeze out distribution functions ff{kx,ky) 
over their transversal and longitudinal momenta kx 
and ky, respectively, at t = 9. Or in RFF. LEFT: 
Evaluation with a vacuum nucleon mass Afjv (0) — 939 
MeV. RIGHT: Evaluation with the density and tem- 
perature dependent nucleon mass A4"jv(n, T). The over- 
all norm of the Freeze out distribution function has 
increased by a factor of ~ 4 when density and temper- 
ature dependent nucleon masses were taken into ac- 
count. 



three-fluid model, which is a highly involved tool, 
cf. [b^. Therefore, for the time being it is difhcult 
to say how strong the impact of mesons is. There- 
fore, we were aiming at a description, which allows 
to account for the nucleon mass shift scenario dur- 
ing the freeze out process in a more common way. 

Finally, we remark that in-medium modifica- 
tions have actually to be taken into account al- 
ready before and during the hadronization pro- 
cess. This points to an even stronger impact of 
in-medium modifications on the final particle spec- 
trum than presented. 



V. SUMMARY 

We have investigated a freeze out scenario within 
a finite layer for a massive nucleon gas. Special 
attention has been drawn about how strong the 
impact of the in-medium nucleon mass modifica- 
tion on the thermal freeze out process is. By fo- 
cussing on a purely nucleon gas we have found a 
substantial effect on the thermodynamical quan- 
tities like temperature T, flow velocity u, parti- 
cle density n and energy density e of the inter- 
acting component. All of these thermodynamical 
functions have revealed a faster freeze out com- 
pared to a scenario without in-medium nucleon 
mass shift. These modifications have a sizable 
implication on the freeze out particle distribution 
function, which is a basic observable in heavy-ion 
collision experiments. For small momenta around 
the nucleon mass a strong change of about a fac- 
tor ~ 2 has been found (see Fig. 2J. A contour 
plot of the particle distribution function in the 
transversal-longitudinal momentum plane (kx,ky) 
illustrates this effect (see Fig. (SJ • From these re- 
sults we conclude that in-medium modifications 
of nucleons have a significant consequence on the 
freeze out process. This reasoning is certainly valid 
for heavy-ion collisions, which produce sufhciently 
high nucleon particle densities; in particular for ex- 
periments like Compressed Baryonic Matter (CBM 
collaboration) planned at the GSI facility in Darm- 
stadt/Germany. 

For a more realistic description of heavy ion col- 
lisions one should include in the analysis at least 
the low laying mesons and baryons as well. All 
hadrons suffer in-medium modifications of their 
masses and widths, but there are strong differ- 
ences among them. For instance, while the pion 
mass remains almost unaffected by the hadronic 
medium even at very high temperatures and den- 
sities, this is not the case for nucleons, kaons and 
Delta resonances. Taking into account the pions 
and the in-medium modifications of other hadrons 
in the fireball produced in nucleus-nucleus reac- 
tions could modify our results in the details, but 
not the general statement that in-medium modifi- 
cations have some relevance for the freeze out pro- 
cess. For example, the implementation of the pi- 
ons leads to a stronger temperature dependence of 
the chiral condensate ^53| , which causes a stronger 
down shift of the nucleon mass with increasing 
temperature. Then, our results might even be 
more pronounced. In addition, the implementa- 
tion of in-medium modifications has to be taken 
into account before and during the hadronization, 
which leads also to an amplification of their impact 
on the whole freeze out process. 

In summary, our findings for a purely nucleon 
gas suggest that taking into account in-medium 
modifications of nucleons seems to be a necessary 
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and interesting phenomenon, in particular for col- 
lision scenarios with high baryonic densities. 
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Appendix 

The function {n = 1,2) are defined as 



G^{Mn,v,T) 



^ jdkk[^k^+Ml 



xEi( ^^/P + A^i^"'^ 



T 



(36) 

where Ei is a special case of incomplete Gamma- 
function l6a| and defined as 



Ei(a;) 



dt 



(37) 



The function Kn is the Bessel function of second 
kind ,66.] . defined as 



2"n! 
(2n)! 



Finally, we prove the vanishing of the second term 
in Eq. Q. First, we note explicitly the relevant 
four-current and energy- momentum tensor compo- 
nents as deduced directly from the microscopic ki- 
netic definitions |2Jl and respectively. We re- 
call that due to the immediate re-thermalization 



limit during the freeze out there is actually a 
Jiittner type distribution for /i, but with the 
thermodynamical functions T and v as evaluated 
with the approach presented and given in Figs. [3 
Therefore, at the beginning of the time-step for 
fi one has to insert the Jiittner distribution H15() . 
but with the evaluated functions T, and v, into 
the microscopic definitions (|2Jl and ||2Jl, getting the 
following components in RFF: 

7V° = ^[2ahKQ{a)+AbKi{a)] , (39) 



= -[2vabKa{a) +4vbKi{a)] , (40) 
T°° = '^[2ab^Ki{a) + 2b^{3 + v^)K2{a)] , 

(41) 



+ - 



nT 



2a b^ V Ki{a) + 2— (1 + 3v^)K2{a) 



(42) 



T 

-2 — N' 
nT 



2ab^v^ Ki{a) + 2b\Z + v^)K2{a) 



(43) 

We recall that a — M/T, and b — -fa with 7 = 
(1 — w^)^^/^. The second term in Eq. is given 
as 

dUf,T'"'u^ = duaT°°Uo + duoT°^u^ 

+ du^T'^^uo + du^ T^'^'u^ . (44) 

With = 7(1, — 0, 0) we get d-uo — "f'^ vdv and 
dux — —7'^ dv. By using these relations and in- 
serting the components (|39|l - H43|l into 144|l we 
immediately find du^T^^^u^ — 0. We recall that 
aK2{a) = aKo{a) + 2Ki{a). 



[1] E. Shuryak, Phys. Rev. Lett. 68 3270 (1992). 

[2] T.S. Biro, E. van Doom, B. Miiller, M.H. Thoma, 

X.-N. Wang Phys. Rev. C 48 1275 (1993). 
[3] X.N. Wang, Nucl. Phys. A 590 47c (1995). 
[4] K. Geiger, J.I. Kapusta, Phys. Rev. D 47 4905 

(1993). 

[5] P. Koch, B. MiiUer and J. Rafelski, Phys. Rep. 

142 167 (1986). 
[6] K. Geiger, Phys. Rep. 258 237 (1995). 
[7] J. Cleymans, K. Redhch, Phys. Rev. C 60 054908 



(1999). 

[8] F. Becattini, M. Gazdzicki, A. Keranen, J. Man- 
ninen, R. Stock, Phys. Rev. C 69 024905 (2004). 
[9] F. Becattini, J. Cleymans, A. Keranen, E. Suho- 
nen, K. Redlich, Phys. Rev. C 64 024901 (2001). 
[10] F. Becattini, M. Gazdzicki, J. Sollfrank, Eur. 

Phys. J. C 5 143 (1998). 
[11] J. Cleymans, H. Satz, Z. Phys. C 57 135 (1993). 
[12] P. Braun-Munzinger, D. Magestro, K. Redlich, J. 
Stachel, Phys. Lett. B 518 41(2001). 



12 



[13] P. Braun-Munzinger, I. Heppe, J. Stachel, Phys. [42] 

Lett. B 465 15 (1999). 
[14] J. Letessier, J. Rafelski, nucl-th/0504028 (2005); [43 

W. Florkowski, W. Broniowski, M. Michalec, Acta [44' 

Phys. Pol. B 33 761 (2002); 

G. D. Yen, M.I. Gorenstein, Phys. Rev. C 59 2788 [45 
(1999); 

J. Sollfrank, M. Gazdzicki, U.W. Heinz, J. Rafel- [46 
ski, Z. Phys. C 61 659 (1994). 

S.A. Bass, NucL Phys. A 698 164 (2002). [47 
X.-N. Wang, NucL Phys. A 698 296 (2002). 
P. Huovinen, Acta Phys. Pol. B 33 1635 (2002). 
W. Florkowski, W. Broniowski, Phys. Lett. B 477 [48' 
73 (2000). [49 

D. Zschiesche, S. Schramm, J. Schaffner-Bielich, 

H. Stoecker, W. Greiner, Phys. Lett. B 547 7 
(2002). [50' 
T. Yamazaki et. al., Phys. Lett. B 418 246 (1998). 
V. Thorsson, A. Wirzba, Nucl. Phys. A 589 633 [51 
(1995). 

V.L. Eletsky, B.L. loffe, Phys. Rev. Lett. 78 1010 [52 
(1997). 

S. Mallik, S. Sarkar, Phys. Rev. C 69 015204 [53 
(2004). 

F. Laue et al. (for the KaoS collaboration), Eur. 
Phys. J. A 79 397 (2000); Phys. Rev. Lett. 82 
1640 (1999). [54 
R. Barth et al. (for the KaoS collaboration) , Phys. 
Rev. Lett. 78 4007 (1997). [55 
C.L. Korpa, M.F.M. Lutz. Acta Phys. Hung. A 22 
21 (2005). [56 
L Tanihata, S. Nagamiya, S. Schnetzer, H. Steiner, 
Phys. Lett. B 100 121 (1981). [57 
H.W. Barz, L.P. Csernai, W. Greiner, Phys. Rev. 
26 740 (1982). 

Cs. Anderlik, Z.L Lazar, V.K. Magas, L.P. Cser- 
nai, H. Stocker, W. Greiner, Phys. Rev. C 59 388 
(1999). 

V.K. Magas, Cs. Anderlik, L.P. Csernai, F. Grassi, [58' 
W. Greiner, Y. Hama, T. Kodama, Zs. Lazar, H. [59 
Stocker, Phys. Lett. B 459 33 (1999). [60 
C. Anderlik, L. P. Csernai, F. Grassi, et al., Phys. 
Rev. C 59 3309 (1999). 

■V.K. Magas, C. Anderlik, L.P. Csernai et al., [61 
Heavy Ion Phys. 9 193 (1999). 

V.K. Magas, A. Anderlik, Cs. Anderlik, L.P. Cser- [62 
nai, Eur. Phys. J. C 30 255 (2003). 

E. Molnar, L.P. Csernai, V.K. Magas, A. Nyiri, K. [63 
Tamosiunas, arXiv: nucl-th/0503047| 

E. Molnar, L.P. Csernai, V.K. Maga s, Zs.I. Lazar, 
A. Nyiri, K. Tamosiunas, arXiv: nucLth/0503048 
L.P. Csernai, "Introduction to Relativistic Heavy 
Ion Collisions" , John Wiley & Sons, Chichester [64 
(1994) ISBN 471 93420 8. 
J. Kapusta, Phys. Rev. C 29 1735 (1984). 
C. Grant, J. Kapusta, Phys. Rev. C 32 663 (1985). [65 
T.S. Olson, W.A. Hiscock, Phys. Rev. C 39 1818 
(1989). 

F. Jiittner, Ann. Phys. und Chem. 34 856 (1911). [66 
C. Rebbi, "Lattice Gauge Theories and Monte 
Carlo Simulations", Addison Wesley, Reading, 
MA, (1988). 



M.A. Shifman, A.I. Vainshtein, V.I. Zakharov, 
Nucl. Phys. B 147 385 (1979); ibid 448; ibid 519. 

B. L. loffe, Nucl. Phys. B 188 317 (1981). 

E.G. Drukarev, Prog. Part. Nucl. Phys. 50 659 
(2003). 

R.J. Furnstahl, D.K. Griegel, T.D. Cohen, Phys. 
Rev. C 46 1507 (1992). 

T.D. Cohen, R.J. Furnstahl, D.K. Griegel, Prog. 

Part. Nucl. Phys. 35 221 (1995). 

L.S. Celenza, CM. Shakin, "Relativistic Nuclear 

Physics: Theories of Structures and Scattering" , 

World Scientific, Signapore, (1986). 

J.D. Walecka, Ann. Phys. 83 491 (1974). 

A.M. Rakhimov, M.M. Musakhanov, F.C. 

Khanna, U.T. Yakhshiev, Phys. Rev. C 58 1738 

(1998). 

R.J. Furnstahl, C.E. Price, G.E. Walker, Phys. 
Rev. C 36 2590 (1987). 

M. Jaminon, C. Mahaux, Phys. Rev. C 40 354 
(1989). 

M. Lutz, S. Klimt, W. Weise, Nucl. Phys. A 542 
521 (1992). 

S. Zschocke, B. Kampfer, O.P. Pavlenko, Gy. 
Wolf, XL International Winter Meeting on Nu- 
clear Physics, Bormio (Italy), January 20 - 26, 
2002, arXiv: nucLth/0202066 
X. Jin, T.D. Cohen, R.J. Furnstahl, D.K. Griegel, 
Phys. Rev. C 47 2882 (1993). 
T.D. Cohen, R.J. Furnstahl, D.K. Griegel, Phys. 
Rev. C 45 1881 (1992). 

S. Eidelman et al. (Particle Data Group), Phys. 
Lett. B 592 1 (2004) (URL: http://pdg.lbl.gov). 
When a immerical value for the chiral conden- 
sate is fixed, then the current quark mass ruq = 
{mu + md)/2 is determined by the Gell-Mann- 
Oakes-Renner relation: mj/^ — — 2m5(0|g(jr|0); 
/,r = 92 MeV is the pion decay constant and 
m,r = 138 MeV is the pion mass. 

C. Runge, Math. Ann. 46 167 (1895). 

M.W. Kutta, Z. L Math. u. Phys. 46 435 (1901). 
J.D. Lambert, D. Lambert, "Numerical Methods 
for Ordinary Differential Systems: The Initial 
Value Problem", New York, Wiley (1991). 
P. Braun-Munzinger, J. Stachel, J. P. Wessels, N. 
Xu, Phys.Lett. B 344 43 (1995). 
J. Cleymans, H. Oeschler, K. Redlich, Phys. Rev. 
C 59 1663 (1999). 

At point A of the freeze out front the corre- 
sponding values are da^^^ = (1,0,0,0), fc^pp = 
(i?fe, |k| cosOfc,0, 0), where Qk is the angle be- 
tween the Freeze out normal three vector da and 
the three-momentum k of the particle. 
E.L. Bratkovskaya, W. Gassing, C. Greiner, M. 
Effenberger, U. Mosel, A. Sibirtsev, Nucl. Phys. 
A 681 84 (2001). 

L. P. Csernai, I. Lovas, J. A. Maruhn, A. Rosen- 
hauer, J. Zimanyi, W. Greiner, Phys. Rev. C 26 
149 (1982). 

M. Abramowitz, LA. Stegun, "Handbook of Math- 
ematical Functions" , Dover Publ. Inc., New York, 
1970, ISBN 486 61272 4. 



